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When comparing recent experimental data from the HI and ZEUS Collaborations 
at HERA for diffractive dijet production in deep-inelastic scattering (DIS) and 
photoproduction with next-to-leading order (NLO) QCD predictions using diffrac¬ 
tive parton densities, good agreement is found for DIS. However, the dijet photo¬ 
production data are overestimated by the NLO theory, showing that factorization 
breaking occurs at this order. While this is expected theoretically for resolved pho¬ 
toproduction, the fact that the data are better described by a global suppression 
of direct and resolved contribution by about a factor of two comes as a surprise. 
We therefore discuss in some detail the factorization scheme and scale dependence 
between direct and resolved contributions and propose a new factorization scheme 
for diffractive dijet photoproduction. 


1. Diffractive ep Scattering 

It is well known that in high-energy deep-inelastic ep-collisions a large frac¬ 
tion of the observed events are diffractive. These events are defined experi¬ 
mentally by the presence of a forward-going system Y with four-momentum 
Py, low mass My (in most cases a single proton and/or low-lying nucleon 
resonances), small momentum transfer squared t = (p — py)^, and small 
longitudinal momentum transfer fraction xp = q{p — PY)/qp from the in¬ 
coming proton with four-momentum p to the system X (see Fig. 1). The 
presence of a hard scale, as for example the photon virtuality = —q^ in 
deep-inelastic scattering (DIS) or the large transverse jet momentum p^p in 
the photon-proton centre-of-momentum frame, should then allow for cal¬ 
culations of the production cross section for the central system X with the 
known methods of perturbative QCD. Under this assumption, the cross sec¬ 
tion for the inclusive production of two jets, e+p —> e-|-2 jets-l-V'-t-F, can be 
predicted from the well-known formulae for jet production in non-diffractive 
ep collisions, where in the convolution of the partonic cross section with the 
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Figure 1. Diffractive scattering process ep eXY, where the hadronic systems X and 
Y are separated by the largest rapidity gap in the final state. 


parton distribution functions (PDFs) of the proton the latter ones are re¬ 
placed by the diffractive PDFs. In the simplest approximation, they are 
described by the exchange of a single, factorizable pomeron/Regge-pole. 

2. Diffractive Parton Distribution Functions 

The diffractive PDFs have been determined by the HI Collaboration at 
HERA from high-precision inclusive measurements of the DIS process ep —> 
eXY using the usual DGLAP evolution equations in leading order (LO) and 
next-to-leading order (NLO) and the well-known formula for the inclusive 
cross section as a convolution of the inclusive parton-level cross section with 
the diffractive PDFs A similar analysis of inclusive measurements has 
been published by the ZEUS Collaboration A longer discussion of the 
extraction of diffractive PDFs can be found elsewhere 

3. QCD Factorization in Hard Diffraction 

For inclusive diffractive DIS it has been proven by Collins that the formula 
referred to above is applicable without additional corrections and that the 
inclusive jet production cross section for large can be calculated in terms 
of the same diffractive PDFs The proof of this factorization formula, 
usually referred to as the validity of QCD factorization in hard diffraction, 
may be expected to hold for the direct part of photoproduction ~ 0) 
or low-Q^ electroproduction of jets However, factorization does not hold 
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for hard processes in diffractive hadron-hadron scattering. The problem is 
that soft interactions between the ingoing two hadrons and their remnants 
occur in both the initial and final state. This agrees with experimental mea¬ 
surements at the Tevatron Predictions of diffractive dijet cross sections 
for pp collisions as measured by CDF using the same PDFs as determined 
by HI ^ overestimate the measured cross section by up to an order of mag¬ 
nitude This suppression of the CDF cross section can be explained by 
considering the rescattering of the two incoming hadron beams which, by 
creating additional hadrons, destroy the rapidity gap ®. 

4. Factorization Breaking in Diffractive Photoproduction 

Processes with real photons ^ 0) or virtual photons with fixed, but low 
involve direct interactions of the photon with quarks from the proton 
as well as resolved photon contributions, leading to parton-parton interac¬ 
tions and an additional remnant jet coming from the photon (for a review 
see ®). As already said, factorization should be valid for direct interactions 
as in the case of DIS, whereas it is expected to fail for the resolved pro¬ 
cess similar as in the hadron-hadron scattering process. In a two-channel 
eikonal model similar to the one used to calculate the suppression factor 
in hadron-hadron processes ®, introducing vector-meson dominated photon 
fluctuations, a suppression by about a factor of three for resolved photopro¬ 
duction at HERA is predicted Such a suppression factor has recently 
been applied to diffractive dijet photoproduction and compared to 
preliminary data from HI and ZEUS While at LO no suppression 
of the resolved contribution seemed to be necessary, the NLO corrections 
increase the cross section significantly, showing that factorization breaking 
occurs at this order at least for resolved photoproduction and that a sup¬ 
pression factor R must be applied to give a reasonable description of the 
experimental data. 

5. Factorization Scale Dependence for Real Photons 

As already mentioned elsewhere ^142^ describing the factorization break¬ 
ing in hard photoproduction as well as in electroproduction at very low 
by suppressing the resolved contribution only may be problematic. 
An indication for this is the fact that the separation between the direct 
and the resolved process is uniquely defined only in LO. In NLO these two 
processes are related. The separation depends on the factorization scheme 
and the factorization scale M^. The sum of both cross sections is the only 
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physically relevant cross section, which is approximately independent of the 
factorization scheme and scale As demonstrated in Fig. 2, multiplying 



Figure 2. Photon factorization scale dependence of resolved (dashed) and direct (dot¬ 
ted) contributions to the diffractive dijet photoproduction cross section (full curve). Also 
shown is the sum of the direct and suppressed resolved contribution (dot-dashed curve). 


the resolved cross section with the suppression factor R = 0.34 destroys 
the correlation of the M-^-dependence between the direct and resolved part 
and the sum of both parts has a stronger M^-dependence than for the 
unsuppressed case {R = 1), where the M-^-dependence of the NLO direct 
cross section is compensated to a high degree by the M-^-dependence of the 
LO resolved part. 

The introduction of the resolved cross section is dictated by perturba¬ 
tion theory. At NLO, collinear singularities arise from the photon initial 
state, which are absorbed at the factorization scale into the photon PDFs. 
This way the photon PDFs become M^-dependent. The equivalent M^- 
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dependence, just with the opposite sign, is left in the NLO corrections to 
the direct contribution. With this knowledge, it is obvious that we can 
obtain a physical cross section at NLO, i.e. the superposition of the NLO 
direct and LO resolved cross section, with a suppression factor R < 1 and 
no M.^-dependence left, if we also multiply the In M.^-dependent term of 
the NLO correction to the direct contribution with the same suppression 
factor as the resolved cross section. We are thus led to the theoretical con¬ 
clusion that, contrary to what one may expect, not all parts of the direct 
contribution factorize. Instead, the initial state singular part appearing 
beyond LO breaks factorization even in direct photoproduction, presum¬ 
ably through soft gluon attachments between the proton and the collinear 
quark-antiquark pair emerging from the photon splitting. This would be in 
agreement with the non-cancellation of initial state singularities in diffrac¬ 
tive hadron-hadron scattering 


6. The Transition Region of Virtual Photoproduction 

We now present the special form of the InAf^-term in the NLO direct 
contribution and demonstrate that the AL^-dependence of the physical cross 
section cancels to a large extent in the same way as in the unsuppressed case 
{R = 1). These studies can be done for photoproduction cs 0) as well 
as for electroproduction with fixed, small Q^. Since in electroproduction 
the initial-state singularity in the limit ^ 0 is more directly apparent 
than for the photoproduction case, we shall consider in this contribution 
the low-Q^ electroproduction case just for demonstration. This diffractive 
dijet cross section has been calculated recently 

A consistent factorization scheme for low-Q^ virtual photoproduction 
has been defined and the full (direct and resolved) NLO corrections for 
inclusive dijet production have been calculated in In this work we 
adapt this inclusive NLO calculational framework to diffractive dijet pro¬ 
duction at low-Q^ in the same way as in except that we multiply the 
In ATy-dependent terms as well as the resolved contributions with the same 
suppression factor R = 0.34, as an example, as in our earlier work 1142 .is 
The exact value of this suppression factor may change in the future, when 
better data for photoproduction and low-Q^ electroproduction have been 
analyzed. We present the In AL..y-dependence of the partly suppressed NLO 
direct and the fully suppressed NLO resolved cross section da/dQ^ and 
their sum for the lowest bin. 
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The NLO corrections for virtual jet photoproduction have been im- 
plemented in the NLO Monte Carlo program JETVIP and adapted to 
diffractive dijet production in The subtraction term, which is absorbed 
into the PDFs of the virtual photon /a/ 7 ( 2 : 7 , Mj), can be found in The 
main term is proportional to ln(M^/( 3 ^) times the splitting function 

p,.^,w=2iv.o^£i±<yii!. (1) 

where z = P1P2/P0Q G [x] 1] and Qi is the fractional charge of the quark 
Qi- Pi and p 2 are the momenta of the two outgoing jets, and po and q are 
the momenta of the ingoing parton and virtual photon, respectively. Since 
= —q^ M^, the subtraction term is large and is therefore resummed 

by the DGLAP evolution equations for the virtual photon PDFs. After 
this subtraction, the finite term which remains in the matrix 

element for the NLO correction to the direct process has the same M 7 - 
dependence as the subtraction term, i.e. InM-y is multiplied with the same 
factor. As already mentioned, this yields the M 7 -dependence before the 
evolution is turned on. In the usual non-diffractive dijet photoproduction 
these two M 7 -dependences cancel, when the NLO correction to the direct 
part is added to the LO resolved cross section Then it is obvious that 
the approximate M 7 -independence is destroyed, if the resolved cross section 
is multiplied by a suppression factor R to account for the factorization 
breaking in the experimental data. To remedy this deficiency, we propose to 
multiply the In ALy-dependent term in M with the same suppression 

factor as the resolved cross section. This is done in the following way: we 
split M(g 2 )_ into two terms using the scale in such a way that the 
term containing the slicing parameter ps, which was used to separate the 
initial-state singular contribution, remains unsuppressed. In particular, we 
replace the finite term after the subtraction by 


M[Q , R)yia — 



( M^z 

[ptH^-z) 


^ P? ^ 

yZQ'^ + TjsS j 




R 


( 2 ) 


where R is the suppression factor. This expression coincides with the finite 
term after subtraction (see Ref. for i? = 1, as it should, and leaves the 
second term in Eq. (2) unsuppressed. In Eq. (2) we have suppressed in 
addition to ln{M‘^/p^) also the z-dependent term ln(z/(l — z)), which is 
specific to the MS subtraction scheme as defined in The second term in 
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Eq. (2) must be left in its original form, i.e. being unsuppressed, in order 
to achieve the cancellation of the slicing parameter (j/g) dependence of the 
complete NLO correction in the limit of very small or equivalently very 
large s. It is clear that the suppression of this part of the NLO correction 
to the direct cross section will change the full cross section only very little 
as long as we choose ~ The first term in Eq. (2), which has the 
suppression factor R, will be denoted by DIRis in the following. 

To study the left-over M.y-dependence of the physical cross section, we 
have calculated the diffractive dijet cross section with the same kinematic 
constraints as in the HI experiment Jets are defined by the CDF cone 
algorithm with jet radius equal to one and asymmetric cuts for the trans¬ 
verse momenta of the two jets required for infrared stable comparisons with 
the NLO calculations The original HI analysis actually used a symmet¬ 
ric cut of 4 GeV on the transverse momenta of both jets The data have, 
however, been reanalyzed for asymmetric cuts 

For the NLO resolved virtual photon predictions, we have used the PDFs 
SaSlD and transformed them from the DIS.y to the MS scheme as in Ref. 

If not stated otherwise, the renormalization and factorization scales at 
the pomeron and the photon vertex are equal and fixed to = Pt jeti- We 
include four flavors, i.e. n/ = 4 in the formula for as and in the PDFs of 
the pomeron and the photon. With these assumptions we have calculated 
the same cross section as in our previous work First we investigated how 
the cross section dcr/dQ^ depends on the factorization scheme of the PDFs 
for the virtual photon, i.e. dcr/dQ^ is calculated for the choice SaSlD and 
SaSlM. Here dcr/dQ^ is the full cross section (sum of direct and resolved) 
integrated over the momentum and rapidity ranges as in the HI analysis. 
The results, shown in Fig. 2 of Ref. demonstrate that the choice of the 
factorization scheme of the virtual photon PDFs has negligible influence on 
da/dQ^ for all considered Q^. The predictions agree reasonably well with 
the preliminary HI data 

We now turn to the M.y-dependence of the cross section with a suppres¬ 
sion factor for DIRis. To show this dependence for the two suppression 
mechanisms, (i) suppression of the resolved cross section only and (ii) ad¬ 
ditional suppression of the DIRis term as defined in Eq. (2) in the NLO 
correction of the direct cross section, we consider da/dQ^ for the lowest 
Q^-bin, G [4,6] GeV^. In Fig. 3, this cross section is plotted as a func¬ 
tion of ^ = M-^/p^ in the range ^ G [0.25; 4] for the cases (i) (light full 
curve) and (ii) (full curve). We see that the cross section for case (i) has 
an appreciable ^-dependence in the considered ^ range of the order of 40%, 
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Figure 3. Photon factorization scale dependence of resolved and direct contributions 
to dtr/dQ^ together with their weighted sums for (i) suppression of the resolved cross 
section and for (ii) additional suppression of DIRig, using SaSlD virtual photon PDFs. 


which is caused by the suppression of the resolved contribution only. With 
the additional suppression of the DIRis term in the direct NLO correc¬ 
tion, the ^-dependence of dcr/dQ^ is reduced to approximately less than 
20 %, if we compare the maximal and the minimal value of dcr/dQ^ in the 
considered ^ range. The remaining ^-dependence is caused by the NLO cor¬ 
rections to the suppressed resolved cross section and the evolution of the 
virtual photon PDFs. How the compensation of the M.y-dependence be¬ 
tween the suppressed resolved contribution and the suppressed direct NLO 
term works in detail is exhibited by the dotted and dashed-dotted curves in 
Fig. 3. The suppressed resolved term increases and the suppressed direct 
NLO term decreases by approximately the same amount with increasing 
In addition we show also dcr/dQ^ in the DIS theory, i.e. without subtrac¬ 
tion of any InQ^ terms (dashed line). Of course, this cross section must be 
independent of This prediction agrees very well with the experimental 
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point, whereas the result for the subtracted and suppressed theory (full 
curve) lies slightly below. We notice, that for the additional sup¬ 

pression of DIRis has only a small effect. It increases dcr/dQ^ by 5% only 

24,25 


7. Conclusion 

When comparing experimental data from the HI and ZEUS Collabora¬ 
tions at HERA for diffractive dijet production in DIS and photoproduction 
with NLO QCD predictions using diffractive parton densities from HI and 
ZEUS, good agreement is found for DIS assuming the HI diffractive PDFs. 
However, the dijet photoproduction data are overestimated by the NLO 
theory, showing that factorization breaking occurs at this order. While 
this is expected theoretically for resolved photoproduction, the fact that 
the data are better described by a global suppression of direct and resolved 
contribution by about a factor of two comes as a surprise. We have there¬ 
fore discussed in some detail the factorization scheme and scale dependence 
between direct and resolved contributions and proposed a new factorization 
scheme for diffractive dijet photoproduction. 
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